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Abstract 

Particles that are heavy compared to the electroweak scale (M ^ mw), and that are 
charged under electroweak SU{2) gauge interactions display universal properties such as 
a characteristic fine structure in the mass spectrum induced by electroweak symmetry 
breaking, and an approximately universal cross section for scattering on nuclear targets. 
The heavy particle effective theory framework is developed to compute these properties. 
As illustration, the spin independent cross section for low-velocity scattering on a nu- 
cleon is evaluated in the limit M ^ mw, including complete leading-order matching 
onto quark and gluon operators, renormalization analysis, and systematic treatment of 
perturbative and hadronic-input uncertainties. 
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1 Introduction 



Cosmological evidence for dark matter consistent with thermal rehc Weakly Interacting Mas- 
sive Particles (WIMPs) motivates laboratory searches for such particles interacting with nu- 
clear targets. Search strategies and detection potential are highly dependent on the WIMP 
mass, M, and its interaction strength with nuclear matter. We consider here the class of 
models where the WIMP belongs to an electroweak SU{2) multiplet. This study is motivated 
in part by the observation that an exact discrete parity arises in Standard Model extensions 
involving confined fermions coupled to electroweak SU{2) [1]. The parity ensures stability of 
the lightest pseudo-Nambu Goldstone mode, which is the electrically neutral component of a 
Lorentz scalar, electroweak SU{2) isotriplet [2j. 

Regardless of the origin for such an SU{2) multiplet, e.g. whether it is a composite or 
fundamental particle, universal behavior emerges in the limit where the WIMP mass is large 
compared to the electroweak scale, M ^ mw- The emergence of these universal properties, 
and corrections to them, can be systematically analyzed using techniques of heavy particle ef- 
fective theories [3] . We focus on the case of a real scalar transforming as a triplet of electroweak 
SU{2), although the results extend straightforwardly to arbitrary SU{2) representations, and 
to higher spin particles. 

At energy scales large compared to mw, the new particle is described by an effective heavy 
particle SU{2) gauge theory, 

= (l)l{iv D + . . .)(t), , (1) 

where 0^ is a scalar heavy particle field, is the heavy particle velocity, and is the SU(2) 
covariant derivative. The leading interactions are thus universal, and corrections depending 
on the mass, or other characteristics of the dark matter particle, are suppressed by powers 
of M. In this paper we determine the general structure of the heavy scalar effective theory 
through 0{1/M^). As an illustrative application, we compute the universal cross section for 
low- velocity scattering of 5'?7(2)-charged WIMPs on a nucleon in the limit M ^ mw- We 
present a complete leading order matching onto gluonic operators, renormalization analysis, 
and systematic treatment of perturbative and hadronic-input uncertainties. 

The remainder of the paper is structured as follows. In Section |2| we construct the relevant 
heavy particle effective theory at scales fi ^ rriw, and compute the leading Wilson coefficients 
in a simple model. In Section |3] we consider the operator basis and mass corrections in the 
low energy theory after integrating out scales fi ~ mw- In Section |4] we perform explicit 
matching computations at the scale fi ~ mw- In Section [5] we perform renormalization group 
evolution from /i ~ mw down to low scales /i ~ 1 GeV > Aqcd, where QCD matrix elements 
are estimated. Section [6] presents the cross section for low- velocity scattering on a nucleon. 
Section [7] presents a summary and outlook. 

2 Heavy scalar effective theory: electroweak symmetric 
theory 

Consider the effective theory and matching conditions for a scalar particle of mass M, charged 
under electroweak SU{2). With obvious modifications, the construction applies to general 
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gauge groups. We start by investigating the effective theory at scales mw <^ jj, <^ M, with 
unbroken electroweak gauge symmetry. 



2.1 Lagrangian 

We work in terms of an effective heavy scalar field </)i,(a;), in the isospin J representation of 
SU{2). The covariant derivative is D^ = d^- ig2W^t'}, and = i[D^,D^]/g2 = W^^j is 
the associated field strength. We let gi, g2, gs = g denote the Standard Model U{1)y, SU{2)w 
and SU{3)c gauge coupling constants, respectively. A typical heavy particle momentum can 
be decomposed as 

= Mv^" + A;^ , (2) 

where v'^ is a velocity, = 1, and is a residual momentum. The basis of operators involves 
the perpendicular derivative, 

D'^ = D^" - v^v D . (3) 

Through 0{1/M^), the scalar heavy particle effective theory in the one-heavy-particle sector 
takes the form, 

= <P.[^- ■ D - c,— + + g,cn + ^g2CM 



16M3 ' 16M3 ' 16M3 ' "^^""^ 16M3 



16M3 16M3 ^•^ 16M3 



2 , e^'^P'^v^v^TiiW^^W.^j 



where we have employed appropriate field redefinitions to remove possible redundant operators 
involving factors oiv-D acting on Note that the operators with coefficients c'j^^ through 
violate parity and time reversal symmetries]^ For the effective theory describing a fundamental 
heavy scalar particle, we have Ci = C2 = cai = 1 and cd = cm = ca2 = cas = caa, = c'^i = 
Cyi2 = = = at tree level [3]. We find that the low-energy manifestation of relativistic 
invariance ( "reparameterization invariance" OE]) implies the exact relations, 

Ci = C2 = 1 , Cm = CD . (5) 



Section 2^ provides a nontrivial illustration of the latter relation. 

The complete lagrangian including Standard Model particles and interactions can be writ- 
ten 

£ = £0 + £sM + £0,SM • (6) 



^ Additional CPT violating operators at ©(l/A/^) and 0{1/M^) are constrained by reparameterization 
invariance to have vanishing coefficient. 
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Here £sm is the usual Standard Model lagrangian, and by convention we have included interac- 
tions with in C(f,. So far our discussion applies to a general irreducible SU (2) representation 
for the heavy scalar field Specializing to the case of a real scalar field, necessarily with 
integer isospin, the effective theory is invariant undei[^ 

o -Vf, , ^ . (7) 

It is straightforward to verify that all interactions in are invariant under this transformation. 

In the one-heavy-particle sector, the remaining terms involving the Higgs field H, gauge 
fields, and fermions are [H = iT2H*) 

C^,SM - + ■ ■ ■ + CQ — + cx + cnQ 

UR^iv-DuR dn'^iv-DdR QLHdR + h.c. QlHur + h.c. 

+ CDu r Cod r Cud VJ^ T C//„- 



2 (G) G'^°^G'a/3 2 (G) ^a^^^^^"'^M/g 2 {G) i ^^'''"' '^^'^ju 2 (G) / ^^''^'^^"^M^^va^^pa 
+ . . . 



Terms odd under ([T]) have been omitted. Subleading terms containing only if, 0^, and their 
covariant derivatives are represented by the first ellipsis in (|8|. Terms bilinear in lepton fields, 
and terms bilinear in the hypercharge gauge field are also present in £^,sm but have not been 
written explicitly. Repeated indices a = 1..3 and A = 1..8 imply a sum over gauge generators. 
Reparameterization invariance implies 

cq = cx. (9) 
2.2 Sample matching calculation 

As an illustration of the construction and matching conditions for the heavy particle lagrangian 
£<^, consider the case of a fundamental scalar, ignoring scalar self interactions (i.e., 0^ terms). 
For the matching of the terms containing a single gauge field, we consider the full theory result 
for the amputated three point function. 



where q = p' — and F{q^) is a model-dependent form factor. Setting p^ = p'"^ = M^, 
v'^ = (1,0,0,0), the matching conditions for scalar scattering from a/i = Oor/i = z gauge 




^ For a real scalar field, the effective theory is obtained by introducing in the field redefinition (l>{x) — 
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field read 



F(0) - F\0)q' 



CD 



8M2 



+ ... 



{p + p 



-m 1 - 



4M2 



/2 



F'iO)q^ + ... 



{P + P') 



p^ 



4M2 



p/2 q2 

+ Cm 



8M2 



p/2 _ p2 

8M2 



(cd - CAf ) + . . . . 



(11) 



An explicit computation of one-loop gauge boson corrections, employing dimensional regular- 
ization in o? = 4 — 2e dimensions, yields 



F{q^) = 1 + 



(47r)2 M2 



C2{r) 



2 4, M' 

l + -log — 

3eiR 3 



1 3 1 , M 

^ \ log — 

24eiR 4 12 ^ /i 



+ . 



(12) 

The quadratic Casimir coefficients for the isospin-J and adjoint representations of SU(2) are 
C2{J) = J{J + 1) and C2{G) = 2. From (11) and (12), after effective theory subtractions the 
renormalized coefficients Cni^), CM{f^) in the MS renormalization scheme are found to be 



Cz)(/i) = CA/(/i) 



tt2(/i) 

4:71 



-8J(J + 1) + 12 



/32J(J+ 1) 



( 



log 



M 



(13) 



Matching for a general ultraviolet completion model, and for other effective theory coefficients 
proceeds similarly. 

Our focus will be on the limit M ^ rriw, where all nontrivial matching conditions at the 
scale /i ~ M become irrelevant]^ We leave a detailed investigation of the model-dependent 
form factor and subleading 1/M corrections to future work. Presently we proceed to investigate 
the leading order predictions of the effective theory at scales /i <^ M. 



3 Integrating out the weak scale 

For scattering phenomena at ~keV energy scales of interest to dark matter-nucleus scattering 
search experiments, we should examine the appropriate effective theory far below the elec- 
troweak scale. Let us begin by integrating out the degrees of freedom at the scale mw For 
definiteness we treat the top quark mass rrit and the Higgs boson mass rrih as parametrically 
of the same order as mw- In following sections, we will renormalize to lower energy scales, 
integrating out the remaining heavy quark degrees of freedom as we pass the bottom and 
charm quark thresholds. The remaining hadronic matrix elements may then be evaluated in 
Uf = 3 flavor QCD to obtain cross section predictions. 

^ In particular models with multiple mass scales, 1/M prefactors can be replaced by inverse powers of 
a smaller excitation energy. It is also of interest to investigate whether large anomalous dimensions could 
enhance the coefficients of particular subleading operators. 
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3.1 Mass correction from electroweak symmetry breaking 

We may evaluate the heavy scalar self energy to obtain mass corrections, 



W Z 1 



-^m= + + +•••• (14) 

The shift in mass due to electroweak symmetry breaking appears as a nonvanishing value of 
E(p) at v-p = 0. We find at leading order in the 1/M expansion, and first order in perturbation 
theory, 

— -J + sm — — Jq 



6M = a2mw 



2 



(15) 



In particular, with Q = + Y = ioi Y = 0, the mass of each charged state is lifted 
proportional to its squared charge relative to the neutral component, 

M(Q) - M(Q=o) = a2Q^mw sin^ ^ + 0{1/M) ^ (170 MeV)Q^ . (16) 

Subleading corrections can be similarly evaluated in the effective theory. Since no additional 
operators appear at (9(1/M°), the result (16) is model independent]^ 



3.2 Operator basis 

The effective theory after electroweak symmetry breaking will include: the heavy scalar QED 
theory for each of the electric charge eigenstates, with mass determined as in (|T5|)p| the 
Standard Model lagrangian with W, Z, h,t integrated out; and interactions, 

^ = + -^SM + ^(j)o,SM + • • • , (17) 

where the ellipsis denotes terms containing electrically charged heavy scalars. For the electri- 
cally neutral scalar, 

= <P:,q=o[^v ■ d - ^^^—^ + C?(l/O}0.,Q=o • (18) 

Note that enforcing the reality condition ([T]) implies the vanishing of cd (= cm)- 

Interactions with Standard Model fields begin at order 1/m^. We restrict attention to 
quark and gluon operators (neglecting lepton and photon operators) and again focus on the 
neutral 0t,,Q=o component, dropping the Q = subscript in the following. Mixing with charged 
scalars will become relevant at order l/m^r in nuclear scattering computations; similarly, we 
restrict attention to flavor-singlet quark bilinears, since matrix elements of flavor- changing 
bilinears are suppressed by additional weak coupling factors. Finally, we neglect operators 



* The mass splitting (16) appears in limits of particular models, e.g. [DIZllH]- 

^ We define the pole mass to include the contributions induced by electroweak symmetry breaking, as 
opposed to introducing residual mass terms for different charge eigenstates [S]. 
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+ ... 



Figure 1: Matching condition for quark operators. Double lines denote heavy scalars, zigzag 
lines denote W bosons, dashed lines denote Higgs bosons, single lines with arrows denote 
quarks, and the sohd square denotes an effective theory vertex. Diagrams with crossed W 
lines are not displayed. 



with derivatives acting on (f)^ or involving 75, since these lead to spin-dependent interactions 
that are suppressed for low-velocity scattering. The basis of operators is then 



,,SM 



m 



3 r-uri) 
w 



(0)^(0) (2) ^(2V. 
^Iq ^Iq ^ ^Iq UiiUy^lq 



(19) 



where we have chosen QCD operators of definite spin. 



oft = ni 



a 



O 



A \2 



(20) 



Here A^^B"^ = [A^B'^ + A'^B^)/2 denotes symmetrization. We employ dimensional regu- 
larization with d = 4 — 2e the spacetime dimension. We use the background field method 
for gluons in the effective theory thus ignoring gauge-variant operators, and assume that ap- 
propriate field redefinitions are employed to eliminate operators that vanish by leading order 
equations of motion. The matrix elements of the gluonic operators, O2 , are numerically 
large, representing a substantial contribution of gluons to the energy and momentum of the 
nucleon. To account for the leading contributions from both quark and gluon operators, we 



compute the coefficients 0^2^ through 0{as) and c^i^ through C(a°). 



4 Weak scale matching 

The matching conditions for quark operators in the Uf = 5 flavor theory at renormalization 
scale n = Ht rat ^ raw ~ ""^h are obtained from the diagrams in Fig. ([!]): 



X 



hA 



Xt 



xt 



tD 



(2) 



C 



.(2) 
-ID 



c 



\VtD? 



4(1+Xt)3 

xt{?> + Qxt + 2x 



3(1 + xt) 



(21) 



where subscript U denotes m or c and subscript D denotes s or h. Here C = [7r«|(/i()][J(J- 
l)/2], Xh = rrih/mw and Xt = mt/mw- We ignore corrections of order rriq/mw for q - 
u, d, s, c, b, and have used CKM unitarity to simplify the results. 
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Figure 2: Matching condition onto gluon operators. The notation is as in Fig. [TJ 
Matching conditions onto gluon operators are from the diagrams of Fig. ([2]): 



An 

4:71 



1 3 + 4xt + 2x 

+ 



6(l + x02 



32, fit 

- TT log 

9 mw 



M2 + 3xt} 

— 4 — — log 



9(1+3;^) 



mw{l + Xt) 



4(12xf - 364 + 36a;? - 12x2 + 3xt - 2) 



Hxt - 1)3 
A8x^ + 24a;^ - 104x^ - 35x: 



log 



Xt 



20x2 



l+Xt 

13xi + 18' 



8xt(-3 + 7x2; 

9(X2 - 1)3 ■ 



log 2 



9(x2 



[1+Xt) 



(22) 



There is no dependence of c^2^ or €^2' on CKM matrix elements in the limit of vanishing 
d, s, b quark masses. The renormalized coefficients are computed in the MS scheme. We have 
employed Fock-Schwinger {x ■ A = 0) gauge [10] to compute the full-theory amplitudes for 
gluonic operators in Fig. [2j The effective theory subtractions are efficiently performed in 
a scheme with massless light quarks, using dimensional regularization as infrared regulator. 
We have verified that the same results are obtained using finite masses and taking the limit 
rriq/mw 0. Details of this computation will be presented elsewhere. 



.(2) 



5 RG evolution to hadronic scales 

To account for perturbative corrections involving large logarithms, e.g. as{fio)logmt/ fio, we 
employ renormalization group evolution to sum leading logarithms to all orders. 
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5.1 Anomalous dimensions 

The spin S = and spin S = 2 operators mix amongst themselves, with 



d 



dlog fx 

where 7^^^^'' are {uf + 1) x (uf + 1) anomalous dimension matrices. The leading terms are 



(23) 



7 



(0) 



/ 


\ 


as 


/o 




\ 








~ 4:71 












i^/g)') 




1^32 ■■ 


■ 32 


-2/3o/ 



+ 



/ 64 
' 9 



7 



(2) 



47r 



3 * 



v 



64 
9 



64 
" 9 3 



+ . . . 



(24) 



where /3 = dg/dlog/j, ^ -/3oas/47r, 7^ = log m^/rf log yU ^ -8a^/47r, 7^^ = gd-frn/dg, 
{P/g)' = gd{(3/g)/dg, and /^o = H — |^/- It is straightforward to include subleading terms 



for 7(0) Eni [12] and 7(2) 
5.2 Integrating out heavy quarks 

At the scale fi = fib ^ iTib, we match onto an nj = 4 theory containing u,d, s,c quarks. The 
matching equations are 





= cf\f^b) 




log 


mb\ 
f^b) 


a~{0)/ 




= Cfgifib) 


+ 0{d') 










= c^Vfc) 


4a , 
log 

3 ^ 


rrib 
fib 




+ 0(a2) 




= Cfgifib) 


+ 0{d), 









4 , rrib 
l + d{ ll + -log — 
3 /ife 



O(a^) 



(25) 

where q = u,d, s,c and a = as{fib,nf = 5)/47r. Quantities without (with) tilde refer to the 
Uf = 4 [uf = 5) theory. The scheme dependence for heavy quark masses enters at higher order. 

For definiteness we use pole masses for rrib and rric, with values taken from [15j. Following 

( s) 

our power counting scheme, we consider one less order of in the matching for c\g relative 

(s) 

to C2 • For later use in the numerical analysis, we have included NLO QCD corrections in 
the spin-0 matching fi6[ \T7\ . Similar to above, we evolve coefficients in the rif = 4 theory to 
the scale fi = fic ^ "m-c- Finally, we match onto rif = 3 and evolve to a low scale /io ~ 1 GeV 
independent of heavy quark masses. 
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6 Matrix elements and cross section 



Having expressed the lagrangian in terms of operators renormalized at the scale /io ~ 1 GeV, 
we require hadronic matrix elements evaluated at this scale. 

6.1 Hadronic inputs 

Let us define the zero-momentum matrix elements of renormalized operator; 



(7V|0f '^^(/i)|iV) ^ -L (k>^k^ _ aZml) fSW) ' (26) 



where tjin is the nucleon mass. Matrix elements refer to a definite (but arbitrary) spin state 
of the nucleon. 

6.1.1 Spin zero 

We recall that the spin-0 operator matrix elements are not independent, being linked by the 
relation [18] 

m^ = {l- 7„) Y.{N\m,qq\N) + l-^{N\{Gl^f\N) , (27) 



q 



derived from the trace of the QCD energy-momentum tensor. Here N = p or n. Neglecting 
7m, Oi^al) contributions to f3{g), and power corrections in the above formula, the definitions 
(26) ensure that /g!v(/^) ~ 1 — ^q=uds fq^N- Corrections arising from (27) are included in the 



numerical analysis. 

For quark operators, define the scale-independent quantities, 

^■kN = ^ {p\[uu + dd)\p) , So = {p\[uu + dd - 2ss)\p) . (28) 

In the numerical analysis, we will neglect the small contributions proportional to \Vtd\^ and 
|VtsP, so that = cf]. Neglecting also the small contribution [19] {md—mu){p\{uu — dd)\p) ~ 
2 MeV, and using approximate isospin symmetry, we then require, for N = p ot n, 

^iv(/i°i + /i°l)^S.^, m^f^% = (S.^ - So) = . (29) 



We use nonrelativistic normalization for nucleon states, {N{p)\N{p')) — {2-^)^ 6'^ {p — p' ) . 
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Parameter 


Value 


Reference 


\Vtd\ 


~ 


- 


\Vts\ 


~ 


- 


\Vtb\ 


~ 1 


- 


rriu/md 


0.49(13) 


m 


rris/md 


19.5(2.5) 


m 




0.047(9) GeV 


m 




0.050(8) GeV 


m 




0.064(7) GeV 


m 


So 


0.036(7) GeV 


m 


mw 


80.4 GeV 


m 


rrit 


172 GeV 


m 


nib 


4.75 GeV 


m 


rric 


1.4 GeV 


m 


rriN 


0.94 GeV 




asimz) 


0.118 


m 


a2{mz) 


0.0338 


m 



Table 1: Inputs to the numerical analysis. 



We consider "traditional" values E^n = 64 ± 7MeV gSj and Sq = 36 ± 7MeV [24j, but 
investigate also the lattice determinations, E^f^ = 47±9MeV |21] and E'/* = 50±8MeV [22]]^ 
We adopt PDG values [2H| for light-quark mass ratios. A summary of numerical inputs is 
presented in Table [1} 

6.1.2 Spin two 

The matrix elements of spin-two operators can be identified as 

45 (/^) = dx x[q{x, /i) + q{x, /i)] , (30) 
Jo 

where q{x,^) and q{x,fi) are parton distribution functions evaluated at scale /x. Neglecting 



power corrections, the sum of spin two operators in (20) is the traceless part of the QCD 

energy momentum tensor, hence independent of scale we have f(?p{lj) ~ 1 — X]q=uds fofpil^)- 
Using approximate isospin symmetry we set 

J u,n J d,p ' J d,n ^ u,p ' J s,n J s,p " V-^ / 



^The latter quantity arises from a naive averaging of = 31 ± 15MeV [H] and = 59 ± lOMeV 
See also [2ill71[2H]. 
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/i(GeV) 










1.0 


0.404(6) 


0.217(4) 


0.024(3) 


0.36(1) 


1.2 


0.383(6) 


0.208(4) 


0.027(2) 


0.38(1) 


1.4 


0.370(5) 


0.202(4) 


0.030(2) 


0.40(1) 



Table 2: Operator coefficients derived from MSTW PDF analysis [15] at different values of 

Table |2] lists coefficient values for renormalization scales /i = 1 GeV, /i = 1.2 GeV and /i = 
1.4 GeV determined from the parameterization and analysis of [15]. 

6.2 Cross section 

The low- velocity, spin-independent, cross section for WIMP scattering on a nucleus of mass 
number A and charge Z may be written 



TT 



\ZMp + {A - Z)Mn\ 



mlA^ 



TC 



■\Mp\ 



(32) 



where Aip and are the matrix elements for scattering on a proton or neutron respectiveljj^ 
and rrir = Mmj^ / {M + mj^) denotes the reduced mass of the dark- matter nucleus system. As 
described in Section 16.11 Ain ~ -Mp up to corrections from numerically small CKM factors 



and isospin violation in nucleon matrix elements. In the M ^ mj^f limit, the cross section 
scales as A^. At finite velocity, a nuclear form factor modifies this behavior [29] . 

As a numerical benchmark, let us compute the cross section for low-momentum scattering 
on a nucleon for a heavy real scalar in the isospin representation J = 1. Figure |3] displays 
the result, as a function of the unknown Higgs boson mass. Using Table [T| we consider 
separately the "traditional" inputs Sttat and Sq, as well as recent lattice determinations of 
Ejf^ and T}^^. For each case, separate bands represent the uncertainty due to neglected 
perturbative QCD corrections, and due to the hadronic inputs. We estimate the impact of 
higher order perturbative QCD corrections by varying matching scales m^/2 < fif < 2mf, 
< fil < 2m^, m^/2 < fi^ < 2ml, 1-0 GeV < fiQ < 1.4 GeV, adding the errors in 
quadrature. 

The renormalization group running and heavy quark matching for spin-2 operators are 
evaluated at LOj^ For spin-0 operators, we find a large residual uncertainty at LO from 
/io, /ic and fib scale variation. The RG running from /ic to /xq from (24) is thus evaluated 
with NNNLO corrections, including contributions to (3/g through O(a^) and 7^ through 
0{af). Accordingly, the spin-0 gluonic matrix element from (27) is also evaluated at NNNLO, 



Explicitly, Mn = fn^^{N\ (Eg=u,d, 



JO)^(O) 
^Iq '^Iq 



^ Up to power corrections and subleading 0{as) corrections, our evaluation is equivalent to an evaluation 
in either the n/=4orn^ = 5 flavors theories, taking the c- and 6-quark momentum fractions of the proton 
as input. We have verified that these results, with the matrix elements taken from [15] . are within our error 
budget. 



11 




m/j(GeV) 

Figure 3: Cross section for low-velocity scattering on a nucleon for a heavy real scalar in the 
isospin J = 1 representation of SU{2). The dark shaded region represents the la uncertainty 
from perturbative QCD, estimated by varying factorization scales. The light shaded region 
represents the la uncertainty from hadronic inputs. 

including contributions to f3/g through 0{a'^) and 7^ through 0{a^). The residual /io scale 
variation is insignificant compared to other uncertainties. We perform the RG running and 
heavy quark matching from fit to fic at NLO. Hadronic input uncertainties from each source 
in Table [1] and Table |2] are added in quadrature. We have ignored power corrections appearing 
at relative order as(mc)AQQj-,/m^; typical numerical prefactors appearing in the coefficients of 
the corresponding power-suppressed operators JTE\ suggest that these effects are small. 

Due to a partial cancellation between spin-0 and spin-2 matrix elements, the total cross 
section and the fractional error depend sensitively on subleading perturbative corrections and 
on the Higgs mass parameter ruh. We find 

apiruh = 120 GeV) = 0.7±0. itol x lO^^^cm^ , ap{mh = 140 GeV) = 2A±0.2tol x 10"^ , 

(33) 

where the first error is from hadronic inputs, assuming S|f* and Sjj^J^ from Table [l| and the 
second error represents the effect of neglected higher order perturbative QCD corrections. For 
the illustrative value mt = 120 GeV, and as a function of the scalar strange-quark matrix 
element E^, we display the separate contributions of each of the quark and gluon operators in 
Fig.i 

7 Summary 

We have presented the effective theory for heavy, weakly interacting dark matter candidates 
charged under electroweak SU{2). Having determined the general form of the effective la- 
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Figure 4: Breakdown of contributions to the matrix element A4p using the representative 
values TUh = 120 GeV and = 47(9) MeV. The labels u^^\ c/(^), s^^) and c/(^) refer to spin-S 
up, down, strange and gluon operator contributions, respectively. The thickness represents 
the la uncertainty from perturbative QCD. The left-hand vertical band corresponds to the 
lattice value S^j^* = 50(8) MeV and the right-hand vertical band corresponds to the range 
Ss = 366(142) MeV deduced from Sttat and Sq in Table [l] 



grangian (|4| through 1/M^, we demonstrated matching conditions for subleading operators in 
a simple model. Using the effective theory, we demonstrated universality of the mass splitting 
induced by electroweak symmetry breaking, and of the cross section for scattering on nuclear 
matter. Subleading terms in the 1/M expansion can be studied systematically using (|4]). 

Our focus has been on the case of an isotriplet real scalar [Ij. For this case, relic abun- 
dance estimates [8] indicate that M < fewTeV in order to not overdose the universe. This 
mass range, combined with the universal cross section, provides a target for future search 
experiments. 

We have presented a complete matching at first nonvanishing order in ag, and at leading 
order in small ratios mw/M, mf,/mw and Aqcd/^^c- We performed renormalization group 
improvement to sum leading logarithms to all orders. The residual dependence on the high 
matching scale /it ~ ~ mw represents uncertainty due to uncalculated higher-order per- 
turbative corrections. Assuming the hadronic input Sl,^* from Table [l| this scale variation is 
the largest remaining uncertainty on the cross section; its reduction would require higher loop 
order calculations. 

Our high-scale matching results for quark operators ( [2T| and spin-zero gluon operators 
agree with mw/M — )■ results presented by Hisano et al. [30], under the identification 
f^t = f^b = l-J'c, i-e., a one-step matching onto the Uf = 3 theory]^ This approach neglects 

To make the comparison to the scattering amphtude for a heavy Majorana fermion with x — X'^i we 
use X = %/2e-™^-="(/i„ + H^) = ^/2e™"•^(/^5 + H^), where K and are spinor fields with (1 - = 
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large logarithms appearing in coefficient functions. The effective theory analysis provides a 
systematically improvable method to resum large logarithms, and to estimate both perturba- 
tive and hadronic-input uncertainties. Our results for matching onto spin-two gluon operators 
are new. 

With obvious modifications, our results can be extended to higher SU (2) representations 
(or SU{2) singlets), fermionic WIMPs, and models with additional low-energy field content 
beyond the Standard Model. The treatment of QCD corrections presented here can be ap- 
plied to compute scattering amplitudes in related models, e.g. models involving coupling to 
hypercharge, or supersymmetric models [29 | [3T | 132]. 
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Note added. While this paper was in writing, the preprint [33] appeared, mentioning the 
invariance ([T]) for Majorana fermions. 
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